Sign structures for short-range RVB states on small kagome clusters 
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We extended the Schwinger boson construction to obtain wave functions that are resonating 
valence bond (RVB) counterparts of the degenerate coplanar classical states on the Kagome lattice. 
We examined all 84 of them on the 36-site cluster and found that they form a narrow energy band. 
On the 12-site cluster, there are only four such states and their superpositions accurately account for 
the second through fourth exact diagonalization (ED) states, while the ED ground state is accurately 
reproduced by allowing a particular two-vison insertion on top of the q = RVB state. Thus, we 
have established the RVB sign structures for the low energy states on this cluster. 

PACS numbers: 



Introduction. Despite a long, enduring interest in 
the Heisenberg antiferromagnet on the Kagome lattice, 
the exact nature of the ground state has remained 
elusive ji"— Recent numerical studies increasingly point 
to a spin-disordered ground stater'"— and the best vari- 
ational state to this date is realized by the projected 
Dirac spin liquidi^"— In a series expansion study^^ Singh 
et al found a Valence Bond Solid (VBS) state with an 
energy lower than the projected Dirac spin liquid. While 
one cannot directly compare between series expansion 
and variational energies, it nevertheless indicates that the 
true ground state might be close to a VBS instability^^ 
This is confirmed by a density matrix renormalization 
group (DMRG) study which found that a VBS state is 
indeed nearby in a generalized parameter space.— In a re- 
cent paper^i^ Huh et al investigated this possibility by 
considering spinless excitations of Z2 spin liquids, where 
condensation of vortices carrying Z2 gauge flux (visons) 
leads to a VBS phase, while gapped visons correspond 
to a spin liquid. In this work, we pursue a complimen- 
tary perspective by considering projected Schwinger bo- 
son (SB) wave functions with vison insertions. 

From classical to RVB-like states. We construct varia- 
tional wave functions by performing a Gutzwiller projec- 
tion on the ground state of the following Schwinger-boson 
mean field Hamiltoniani^"— : 
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iMAM\ A = diag{Aa}. 

Here the matrix of real- valued "pairing amplitudes" Aij 
is antisymmetric to produce spin-singlet states; the cor- 
responding Hermitian matrix —iA is diagonalized to A 
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1: [Color online] The SB ansatze {Ay } from Ref. Il6l . 
ij — 1 if an arrow points from site i to j, and —1 for the op- 



posite direction. This construction scheme readily generalizes 
to any ABC pattern by pointing arrows A ^ B ^ C ^ A. 



by a unitary matrix M. The Gutzwiller operator Vg en- 
forces the constraint 6^.60- = 1 on each site. In the 
valence bond basis, Eq. can be viewed as a resonating 
valence bond (RVB) wave function with bond amplitudes 
Uij. In the Sz basis used in our variational Monte Carlo 
(VMC), the wave function is given by the permanent of 
a matrix Ui^^-fj^i^ In a recent work^iS we showed that 
the chemical potential fj, determines the range of the RVB 
amplitudes Uij in the resulting spin liquid and performed 
extensive energetics studies on the original SB ansatze 
from Refs. In this paper, we take advantage of the 

flexibility of the VMC setup to explore some questions 
accessible with more general SB constructions, including 
how the classical degeneracy translates to the quantum 
problem and how vison fluctuations may be important. 

We take n and {Aij } as our variational degrees of free- 
dom. For simplicity, from here on we restrict non-zero 
Aij to nearest-neighbor only and taking values ±1. We 
generalize the symmetric q = and -y/S x -y/S ansatz o^^'^^ 
from our recent worki^ to obtain projected wave func- 
tions that do not respect lattice symmetries. This is 
motivated by the observation that the Kagome antifer- 
romagnet has a highly degenerate manifold of classical 
ground states. Among these, coplanar states are believed 



2 



ABC states 



Dirac SL 



ED ground state 



-0.415 



-0.420 



-0.425 



-0.430 



-0.435 



-0.440 



FIG. 2: Variational energies per site of the 84 "ABC" 
Schwinger boson wave functions on the 36-site cluster. We 
include the Dirac SL and ED energies from Ref. [l9ll23 for 
comparison. We set = for each ansatz a in 

the VMC calculations. The narrow band of variational ener- 
gies suggests that a projection of the Hamiltonian onto the 
subspace spanned by these states may allow a significant im- 
provement in the variational ground state energy, but we did 
not pursue this in the present work. 



to be most importanl^i^— and have spins pointing in one 
of the three coplanar directions A, B, and C forming 120° 
angles such that each triangle has precisely one of each la- 
bel. For each classical ground state, we obtain the corre- 
sponding ansatz as follows By drawing arrows between 
nearest-neighbor sites, pointing from A ~^ B ^ C ^ A, 
we define Aij = -1-1 if an arrow points from site i to j, and 
— 1 for the opposite direction. This "ABC" construction 
scheme indeed reproduces the q = and a/3 x \/3 ansatze 
found by Sachdev (see Fig. [T])) ^^-^^ but more generally, 
it also admits spin liquids that break some lattice sym- 
metries and we refer to them as "non-symmetric" . Note 
that permuting the labels A, B, and C does not lead to 
a distinct ansatz since the direction of the arrow linking 
sites i and j remains unchanged under a cyclic permu- 
tation; for non-cyclic permutations, the directions of all 
arrows are reversed, but the resulting ansatz is equivalent 
up to a gauge choice. 

On the fully symmetric 36-site Kagome cluster, there 
are 84 distinct ansatze, and we calculate the nearest- 
neighbor Heisenberg exchange energy of the correspond- 
ing projected Schwinger boson wave functions using 
VMC. Note that for each ansatz, {Ay} are fixed at ±1 
(our simplifying choice), and /i becomes the only remain- 
ing variational parameter. To allow a fair comparison 
between the energies of different ansatze as well as to 
reduce the computational cost, we fix /i„ = 1.05^™^'' 
for each ansatz a, where is the threshold value for 

spinon condensation (it is negative and only n < j/^^^ is 
allowed). Figure [2] shows that their energies lie within 
a narrow band above the exact diagonalization (ED) 
ground state, where the latter energy is taken from past 



B-C-B-C-E 
/ \ / \ / 

AAA 
/ \ / \ / 

B-C-B-C-B 
/ \ / \ / 

AAA 
/ \ / \ / 

3-C-B-C-B 
(a) Ansatz I 



C-A-B-A-C 
/ \ / \ / 
B C B 

/ \ / \ / 

C-A-B-A-C 
/ \ / \ / 

B C B 

/ \ / \ / 

C-A-B-A-C 
(c) Ansatz III 



C-B-C-B-C 
/ \ / \ / 
AAA 
/ \ / \ / 

B-C-B-C-B 
/ \ / \ / 

AAA 
/ \ / \ / 

C-B-C-B-C 
(b) Ansatz II 



A-B-A-C-A 
/ \ / \ / 

B C B 

/ \ / \ / 

A-C-A-B-A 
/ \ / \ / 

C B C 

/ \ / \ / 
A-B-A-C-A 
(d) Ansatz IV 



FIG. 3: [Color online] The four distinct "ABC" ansatze on 
the 12-site lattice. Colored labels mark sites on the cluster. 
Ansatz I can be identified with Sachdev's q = ansatz, while 
Ansatz II is obtained from Ansatz I by shifting the last row 
by half a horizontal lattice vector. The remaining ansatze III 
and IV are related to Ansatz II by 60° rotations about the 
center of a hexagon. 



studiesj^^Ti^i For comparison, we also calculated Dirac SL 
energy for the same 36-site cluster i^°i^^i^^i^^ The vertical 
scale is comparable to Fig. 4 in the ED study by Sindz- 
ingre et al which reveals many ED states populating 
our exhibited window <^ We remark that both the q = 
and V3 X energies can be further lowered by varying 
the chemical potential fi and including further-neighbor 
Aij but the lowest achieved energy within such fam- 
ily is about —0.420 per site and is only a marginal im- 
provement on the scale of Fig. [2] We expect roughly 
similar insensitivity to fi in the other "ABC" ansatze as 
well. 

In a recent study,— we included second-neighbor 
Heisenberg coupling J2 and showed that the q = ansatz 
wins over the Dirac spin liquid for J2 > 0.08. Here we 
also find that the \/3 x \/S ansatz wins over the Dirac SL 
foi' J2 ^ —0.05. The q = and x -s/S states have cor- 
respondingly all antiferromagnetic and all ferromagnetic 
second-neighbor spin correlation, while the other "ABC" 
states straddle the two limits and their energies spread 
out from the narrow band in Fig. [2] upon adding J2 . 

Returning to the most challenging problem with J2 = 
0, one may reason from perturbation physics that non- 
vanishing off-diagonal Hamiltonian matrix elements be- 
tween closely spaced energy levels might allow a substan- 
tial gain in energy. However, projecting the Heisenberg 
Hamiltonian into the subspace spanned by the "ABC" 
Schwinger boson wave functions is costly for the 36-site 
VMC calculations with permanents and is not done here. 

On the 12-site Kagome cluster, there are only four dis- 
tinct ansatze as shown in Figs. |3(a)||3(dy These ansatze 
include the q = ansatz, while the three remaining 
ansatze transform into one another under 60° rotations. 
Due to the small Hilbert space, the Hamiltonian matrix 



3 




(a) 2-vison configuration I (b) 2-vison configuration II 



FIG. 4: [Color online] The solid and open circles respectively 
denote sites in the 12-site cluster and equivalent sites under 
periodic boundary condition. The colored sites emphasize the 
q = translational symmetry of the underlying ABC ansatz. 
Each colored plaquette contains a vison, marked by a cross. 
The Aij on each thick link crossed by the path joining a pair of 
visons is multiplied by —1. Under translations and rotations, 
two-vison configuration I generates 12 distinct states, while 
configuration II generates 24 distinct states. Each set leads 
to a unique linear combination with quantum numbers Reo ~ 
(7 = — 1, Tei = = +1) where a denotes refiection about the 
mirror passing through a hexagon center and a midpoint on 
its edge, and Tg-^ and Tg, are generators of lattice translations. 



elements and the wave function overlaps can be computed 
exactly, and we obtain the eigenvalues for the Heisenberg 
Hamiltonian restricted to this four-dimensional subspace. 

Tabic H] compares the variational energies to the three 
lowest ED energies for the 12-site cluster. The lattice 
symmetries i?6o, tr, T^j^, and T^^ are the same as Ce, cr, 
ai and 02 defined in Fig. 1 of Ref. which refer to a 
60° rotation about a hexagon center, a reflection about 
the mirror passing through a hexagon center and a mid- 
point on its edge, and non-parallel translations that join 
a hexagon center to nearest neighbour hexagon centers. 
The optimized chemical potential Ha ~ 1.01/i™^'' is used 
for each ansatz a. An examination of the quantum num- 
bers shows that the best "ABC" state has captured the 
lowest ED excited state. Our analysis of the quantum 
numbers also shows that no linear combination of the 
four states has an overlap with the ED ground state. To 
see this, we first note that the symmetric ansatz has rota- 
tion quantum number Rqq = -1-1. The rotational proper- 
ties of the remaining ansatze only allow one to construct 
three eigenstates of the rotation operator with eigenval- 
ues -1-1, e^^^'^/s. Thus, the lowest and highest states are 
linear combinations of the two Rqq — +1 states, while 
the remaining states with i?6o = e^'^'^/'^ account for the 
two degenerate trial energies that match reasonably the 
degenerate ED pair in the Table. Since the rotation quan- 
tum number of the ED ground state is —1, there is no 
overlap between the "ABC" subspace and the ground 
state. The "ABC" states therefore realize the ED lowest 
excited states. 

Schwinger boson states with vison fluctuations. Tak- 
ing a cue from the recent work of Huh et al^ we con- 
sider a two-vison insertion on top of the q = ansatz. 
This is realized in our work by flipping the sign of Aij 
on the links cut by a string joining the two visons. To 
construct states with quantum numbers Rqq = ct = — 1, 
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TABLE I: Comparison between energies from ED, projected 
"ABC" wave functions, the q = ansatz with two-vison in- 
sertion, and short-ranged resonating valence bond (SRRVB) 
states from Ref. The degeneracy (Deg) and the quan- 
tum numbers of the states are also tabulated. For the "ABC" 
states with no visons, we used optimized ~ 1.01^™^". Our 
"ABC" energy —0.44397 is much closer to the ED first excited 
level compared to —0.43764 obtained in Ref. [23, thus suggest- 
ing that the system wants longer-ranged bonds which are not 
realized in their nearest-neighbor dimer covers. For the two- 
vison states, the energies are not sensitive to and we set 
fia = —00 for these cases in order to obtain better compari- 
son with the SRRVB states; the first line corresponds to the 
superpositions of Fig. [IJb) and the second line to Fig. [4ja). 



Te^ = = -1-1 on the 12-site cluster, it is sufficient to 
consider the ansatze generated by translations and ro- 
tations on the two-vison configurations shown in Fig. |4l 
All other two-vison configurations do not lead to quan- 
tum numbers of the ED ground state. Under transla- 
tions and rotations, configuration I generates 12 distinct 
states with visons residing on corner-sharing triangles, 
while configuration II gives 24 distinct states with visons 
residing on a hexagon and a second-neighbor triangle (to 
see the distinction in some cases one also needs to con- 
sider loop products of Aij going around the torus). For 
each set of ansatze, we obtain a unique (equal weight) lin- 
ear combination with the same quantum numbers as the 
ED ground state. Thus, we have constructed trial states 
that contain two-vison quantum fluctuations, which ap- 
parently helps to lower the energy of the system. 

Table U shows the energies of the projected Schwinger 
boson wave functions with two-vison insertions on top 
of the q = ansatz on the 12-site cluster. As with the 
"ABC" calculations, {^y} arc fixed at ±1 and there- 
fore fi is the only variational parameter. We further set 
/i — 00 to realize short-ranged RVB statesi^ (the re- 
sults are essentially insensitive to n). Our second vison 
configuration, Fig. Hfb), produces a variational energy 
that is very close to the ED ground state. The over- 
lap with the ED ground state is 99.91%, which is very 
large given that there is no variational parameter. In 
the Table, we also list the energies found in an earlier 
study by Mambrini et al^^ where the authors projected 
the Heisenberg Hamiltonian onto the subspace spanned 
by all possible short-ranged dimcr-covers of the 12-site 
Kagome cluster. In contrast to their variational ground 
state which was obtained by diagonalizing the projected 
Hamiltonian in the large subspace, we emphasize that 
our construction is essentially parameter-free. Further- 
more, it is difficult to understand the physics behind the 
sign structure of the variational ground state in the RVB 
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basis. The simplicity of our trial zero-momentuni state 
and its ability to reproduce their energy suggest a re- 
markable picture of vison dynamics hidden within the 
complex picture of randomly fluctuating valence bonds. 

Discussion. The strong overlap between our zero- 
parameter projected Schwinger boson wave function with 
the two-vison insertion and the ED ground state on the 
12-site cluster adds some support to the increasingly com- 
mon view that the ground state of the Kagome antifer- 
romagnet is close to a VBS transitioni^ i^^i^^ We can con- 
template a picture where the ground state contains sig- 
nificant local two-vison fluctuations (note that as long as 
individual visons are not free, the system remains a spin 
liquid). It is difficult to construct workable wave func- 
tions with such local fluctuations on larger clusters, since 
we would expect the number of such desired pair-vison 
insertions to grow with the system size. If we take the 12- 
site study with one pair-vison as a guide, on the 36-sitc 
cluster we would want on the order of three pair-vison in- 
sertions and correspondingly complex superpositions to 
produce a symmetric trial state. 

So far on the 36-site cluster, we have only considered 
single pair-vison insertions and constructed superposi- 
tions that produce translationally invariant states. For 
a particular local pair-vison picture, we sum over lattice 
translations and rotations of the corresponding SB wave 
functions to obtain the trial state. Such superpositions of 
permanents are computationally costly, but we can still 
report few results. 

We again take the q = SB state with no visons as 
the starting point. Adding a static two-vison defect leads 
to the increase in energy; for example, for wave functions 
with two-vison insertions shown in Fig. 3] the trial energy 
became -0.403(2) per site (an increase by about 0.015 per 
site or 0.5 in total energy from the no-vison case, which 
gives some idea about the vison core energy). When we 
tried translationally invariant superpositions, we found 
that the energy further increased for Fig. SJa) or re- 
mained unchanged for Fig. HKb), which is qualitatively 
different from the 12-site cluster where we observed sig- 
nificant lowering of energy. We note, however, that the 
tried states have the mirror and rotation quantum num- 
bers a = —I, Rea = —1, which are different from the 
quantum numbers of the true ground state on the 36-site 
clusterj^ Thus, the defect constructions that were opti- 



mal on the 12-sitc cluster are not so on the 36-site cluster, 
indicating that finite-size effects are still important in the 
energetics and that we should not take the small cluster 
findings too literally. The increase in energy of the zero- 
momentum superpositions of these defects suggests that 
the specific pair-vison dispersion has a band minimum 
at a different momentum, which we would like to exam- 
ine in the future. We have also tried several pair-vison 
states that have the mirror and rotation quantum num- 
bers a = -|-l,i?60 = +1. For example, in the case where 
the visons are located on a hexagon and a triangle that 
share a bond, the trial energy of the wave function with 
static defect is -0.412 per site and decreased to -0.416 
per site upon forming a translationally invariant super- 
position. The latter value is close to the best SB wave 
function energy with no defects. It would be interest- 
ing to explore other superpositions and also think about 
states with several pair-vison fiuctuations. 

This describes our limited work so far, but it shows 
potential for more detailed studies of the properties of 
the SB spin liquid and its excitations. For example, we 
can contemplate estimating vison hopping amplitudes by 
calculating transition matrix elements for moving a single 
vison, which arc important ingredients in the analysis of 
possible vison condensation patterns in Ref. [isl 

We conclude by highlighting possible extensions. As 
already mentioned, it would be interesting to project the 
Hamiltonian into the space spanned by the 84 "ABC" 
states on the 36-site cluster. The needed matrix el- 
ement calculations are difficult in VMC but probably 
within reach, particularly if utilizing the expected lat- 
tice symmetries. Using superpositions of permanents, 
one could also study SB wave functions on clusters with 
odd number of sites like the symmetric 27-site Kagome 
cluster^ It would also be interesting to measure energies 
for lattice-symmetric superpositions in the slave-fermion 
approaches, for example for the popular Dirac and neigh- 
boring Z2 spin liquids ^^—1^ and also allow two-vison in- 
sertions like we did for the SB construction here. Finally, 
while we focused on the J2 = model, it would be in- 
teresting to follow the energetics upon adding small J2 
where we expect the SB states to perform wellJ^'^^i^'^ 
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